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Abstract 

A combined fit is performed to the BaBar and Belle measurements of the e~^e~ tt~^tt~^{2S) 
cross sections for center-of-mass energy between threshold and 5.5 GeV. The resonant parameters 
of the y(4360) and y(4660) are determined. The mass is 4355t?o ^ 9 MeV/c^ and the width 
is 1031^5 =t 11 MeV/c^ for the y(4360), and the mass is 466lt8 ± 6 MeV/c^ and the width is 
42+^^ ± 6 MeV/c2 for the y(4660). The production of the y(4260) in -K+Tr-^ilS) mode is found 
to be at 2(7 level, and ^(y(4260) tt+tt-'iI){2S)) ■ Tg+g- is found to be less than 4.3 eV/c^ at 
the 90% confidence level, or equal to 7.41^;^ eV/c^ depending on it interferes with the y(4360) 
constructively or destructively. These information will shed light on the understanding of the 
nature of the Y states observed in initial state radiation processes. 
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I. INTRODUCTION 



Charmonium spectroscopy is of great interest to both the experimentahsts and the the- 
orists since its first discovery more than 30 years ago. With the successful running of the 
two S-factories at SLAC and KEK, charmonium physics was also revitalized as more and 
more charmonium and charmonium-like states were observed in B decays, in initial state 
radiation (ISR) processes, in double-charmonium productions, and in two-photon processes. 
The observation of the X(3872) Q, y(4260) 0,11], Z(4430)+ 8 and so on may suggest the 
existence of new type of hadronic states besides the conventional mesons (qq) and baryons 
(qqq) in quark model. Among these many newly observed states, those with J^^ = 1 are 
becoming very puzzling since there are too many of them between 4-5 GeV/c^ than expected 
from the potential models J5|,|6|. In addition to the excited ip states observed in the inclusive 
hadronic cross section (the V'(4040),i;(4160), ^(4415)), there are four new structures 

observed in ISR processes (the F(4008) [sj, F(4260) 0, Q, ^(4360) F(4660) @). 

The overpopulation of the vector states in this mass range may suggest at least one of them 
are non-conventional charmonium state 11 ll. 



In the measurement of the 7r^7r~ip{2S) cross section via ISR at BaBar [10|, a broad reso- 
nance like structure was observed with mass 4324 ±24 MeV/c^ and width 172 ±33 MeV/c^; 
while with more luminosity, the Belle experiment found it is indeed due to two narrow res- 
onances with masses of 4361 MeV/c^ and 4664 MeV/c^ 0]. A close examination of the 
BaBar observation shows that there is one bin with high cross section at around the second 
resonance observed in Belle experiment although it is not very significant. 

In this paper, we try to perform a combined fit to e"^e" tt^'k~iI){2S) cross sections 
measured by the BaBar and Belle experiments with two resonances, to obtain a better 
estimation of the resonant parameters of the y(4360) and y(4660). We also study the 
possible production of the F(4260) in ■n^'K~ip{2S) mode. 



II. THE DATA 



Both BaBar and Belle experiments reported cross sections of e^e~ 7r+7r~?/^(2S') for 
center-of-mass energy ranges from threshold to 5.5 GeV. The integrated luminosity of the 
BaBar data sample is 298 fb~^ while that of the Belle data sample is 673 fb~^, with ~ 90% 
of the data were collected at the T(4S') resonance {^/s = 10.58 GeV), while the rest were 
taken off the T(4S') peak. 

Figure [1] shows the data, good agreement between BaBar and Belle results is observed, 
and the two structures are evident. 



III. THE FORMULAE AND THE LIKELIHOOD FIT 



The Breit-Wigner form of a single resonance used in this analysis is the same as in 
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FIG. 1: The cross sections of e~^e~ — > ^^~^^T~^lJ{2S) measured at BaBar (stars with error bars) 
and Belle (dots with error bars), (a) shows all the data, while (b) shows the region with the Y 
resonances. 



where M is the mass of the resonance, Ftot and Fg+e- are the total width and partial width 
to e~^e~ respectively, B{R f) is the branching fraction of R decays into final state /, and 
PS{^/s) is the three-body decay phase space factor. 

In fitting to the data, we assume all the cross sections are due to the y (4360) and F(4660) 
resonances, and they are added coherently, that is, 

a(v^) = I Wi(v^) + BW2{y^) ■ e'^W (2) 

where BWi and BW2 represent the two resonances and is the relative phase between them. 

We fit the data using a binned maximum likelihood method with MINUIT in the CERN 
Program Library [l2|. To take the Poisson distribution of the small number of events in 
each 7T~^7r~ip{2S) mass bin into consideration, we start from the observed number of events 
in each mass bin instead of from the measured cross section. For each 7T~^7r~tlj{2S) mass 
bin, the probability of observing n°^^ events when the expected number of events is Aj is 

ji = ' „obsi • -tiere 

Xi = aieiCiB{^{2S) ^ n+ n- J / ^)B{J / ^ tT) + n^^^, (3) 
where cxj is the mean cross section in the bin (with bin width Am), which is calculated 
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FIG. 2: The results of the binned maximum hkehhood fit to e'^e~ — > 7r^7r~^/^(25') data from Belle 
and BaBar. The curves show the best fit with two coherent Breit-Wigners and the contribution 
from each component. The interference between the two amplitudes is not shown. The two dashed 
curves at each peak show the two solutions (see text). 
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the number of background events in each bin, n 
Belle, and n^^^ is neglected for BaBar. eiCiB{ip{2S) 
bin is obtained through [n"^^ — n°^^) / a^^^ from the vr 
and cross sections presented in Refs. jl, 13 • 

The likelihood is defined as the product of fi over all the 7r~^7r~ilj{2S) mass bins, that is 
£ = rii fi- In reality, — 21n£ is minimized to get the best estimation of the parameters. 



IV. THE y(4360) AND y(4660) 



We fit the Belle data on 7r^7r~ip{2S) between threshold and 5.5 GeV (in 60 bins) and the 
BaBar data in the same energy range (in 30 bins) simultaneously. Figure [2] shows the fit 
results. Two equally good solutions are found with the two amplitudes interfere with each 
other differently. The masses and the widths of the two resonances are identical but the 
partial widths to e^e~ and relative phases are different in these two solutions, as shown in 
Table [H The statistical significance of the y(4660) is calculated by comparing the likelihood 
of the fit with and without it, and we found it is 6. la in the combined fit, while that quoted 
in Belle experiment is 5.8cr |9|. 

To validate our fitting method, we fit the Belle data only with two coherent Breit-Wigners, 
and the results are shown in Table [III together with the Belle results j9| from an unbinned 
maximum likelihood fit method. From the table, we found that while the two fits agree with 
each other reasonably well, the mass of the second resonance shifts by about 9 MeV/c^, 
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TABLE I: Fit results to the combined BaBar and Belle data on e~^e — > vr+vr '0(25'). The errors 
are statistical only. M, Ttot) and BT^+g- are the mass (in MeV/c^), total width (in MeV/c^), 
and product of the branching fraction to 'K~^TT~ip{2S) and the e~^e~ partial width (in eV/c^), 
respectively, (p is the relative phase (in degrees). 



Parameters Solution I Solution II 
M(y(4360)) 4355+'Jo 



rtot(^(4360)) W3t\l 

i3re+e-(y(4360)) n.itli 12.3 ±1.2 

M(y(4660)) 46611^ 
rtot(l^(4660)) A2t\l 

Br^+^-{Y{AmO)) 2.2toi 5.9 ±1.6 



TABLE II: The binned fit to Belle data. The errors are statistical only. The definitions of the 
parameters and the units are the same as in Table [H Numbers in parentheses are taken from 
Ref. 0] from an unbinned fit. 

Parameters Solution I Solution II 

M(y(4360)) 43591^^0 (4361 ± 9) 

rtot(l^(4360)) 85tJ^(74±15) 

^Pe+e- (1^(4360)) 11.9ti:^ (10.4 ± 1.7) 12.8tJ:^ (11.8 ± 1.8) 

M(y (4660)) 46551^^ (4664 ± 11) 

Ptot(l^(4660)) 40;j^(48±15) 

OTe+e- (y(4660)) 3.4t^;^ (3.0 ± 0.9) 6.0l3;5 (7.6 ± 1.8) 

(j) itjt (39 ± 30) (-79 ± 17) 



which is comparable to the error quoted by Belle. One of the reasons of the shift is the 
fitting method, that is, due to the binning of the experimental data. 

In principle, the difference between a binned and an unbinned likelihood fit should be 
very small when the data sample is large enough. However, in our special case, the data 
sample is rather small (62 events in BaBar experiment and 110 events in Belle), and thus the 
fluctuation due to binning is not small. In other words, the lost of information in the binning 
process can not be neglected. We test this by doing toy experiments with Monte Carlo (MC) 
simulation. We generate 100 MC samples each with the same number of observed events as 
in Belle experiment (110 events) j^. For each sample, we do binned and unbinned likelihood 
fits and check the difference between the fit values of the parameters. Figure [3] shows the 
distributions of the differences in masses and widths from the two kinds of fit, and Table UTTl 
shows the fit results to the distributions with Gaussian functions. 

It can be seen that the mean values of the differences are consistent with zero as ex- 
pected, while the standard deviations are significantly different from zero, indicating that 
the expected uncertainties introduced by using a binned fit are at a few MeV/c^ level from 
the unbinned fit. The observed difference in the mass of the y(4660) between binned and 
unbinned fits can be explained as the uncertainty introduced by the binning procedure. We 
take this kind of differences as one extra source of the systematic error. It should be noted 
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FIG. 3: The distributions of the parameter difference between binned and unbinned hke- 
hhood fits. The curves show the fits with Gaussian functions. (a) M^™(y(4360)) — 
M^'*"^(y (4360)), (b) r|?^^(y(4360)) - r^„f ''^(y(4360)), (c) Afb"^(y(4660)) - M"°bi°(y (4660)), 
and (d) rbin(y(4660)) -r-fH^(4660)). 

TABLE III: Expected differences between binned and unbinned fits. AM and AFtot are the 
differences in mass (in MeV/c^) and in total width (in MeV/c^). /i is the mean value and a is the 
standard deviation of a Gaussian distribution, the errors are statistical. 



Parameters 






a 




AM(y(4360)) 


0.3 ±0.3 


2.6 


± 


0.3 


AFtot(^(4360)) 


0.2 ±0.6 


4.4 


± 


0.5 


AM(y(4660)) 


-0.8 ±0.4 


3.0 


± 


0.5 


AFtot(>^(4660)) 


-1.0 ±0.7 


5.2 


± 


0.6 



that the correlations between the parameters are neglected and we do not try to quantize the 
uncertainties of the B ■ Fg+e- measurements here to simplify the error estimation procedure. 

The other sources of the systematic errors on the mass and width measurements are 
not very different from those listed by the Belle experiment 0, since it is mainly due to 
the parametrization of the resonance. Taking directly the error from Belle and adding 
in quadrature with the uncertainty due to the binning, one gets the total systematic er- 
rors of 9 MeV/c2, 11 MeV/c2, 6 MeV/c^, and 6 MeV/c^, for M(y(4360)), Ftot(>"(4360)), 
M(y(4660)), and Ftot(>"(4660)), respectively. 
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FIG. 4: The results of the fit to e^e~ — > Tr'^7r~ip{2S) data from Belle and BaBar. The solid curve 
show the best fit with three coherent Breit-Wigners: the y(4260), y(4360), and y(4660), and the 
dashed curve is the signal shape of the y(4260). 

V. THE y(4260) 

The y(4260) was observed in ir^ir^J/ip mode 0, [sj, and its production in 7r~^7r~ip{2S) is 
not forbidden by any selection rule, except that the phase space is a bit small due to large 
iIj{2S) mass. We add the y(4260) amplitude in Eq. [2] coherently with mass (4247 MeV/c^) 
and width (108 MeV/c^) fixed to the Belle measurement to measure its production rate. 

Figure m shows the fit result with three resonances, the change of — 21n£ is 5.2 for two 
more free parameters, this corresponds to a statistical significance of the y(4260) of 1.8a. 
That is, the production of the y(4260) in 7r~^7r~ilj{2S) is not significant. There are two 
pairs of solutions for ,B(y(4260) 7r+7r~?/^(2S')) ■ Fg+g-, one pair interferes with the y(4360) 
constructively, with i3(y(4260) tt^tt" -?^(2^))-Fe+e- = 1.4to:9 eV/c^, and we set the upper 
limit as 4.3 eV/c^ at the 90% confidence level; the other pair interferes with the y(4360) 
destructively, and i3(y(4260) ^ n+7i-ij{2S)) ■ T^+e- = 7.4^?:^ eV/c^. 



VI. SUMMARY 

Based on combined Belle and BaBar data on e+e~ —>■ 7r^7r~ip{2S), we perform a 
maximum likelihood fit to get the resonant parameters of the y(4360) and y(4660). 
We obtain M(y(4360)) = 4355l?o ± 9 MeV/c^, Ftot (1^(4360)) = 103t\l ± H MeV/c^, 
M(y(4660)) = 466li| ± 6 MeV/c^, and Ftot(l'(4660)) = 42t\l ± 6 MeV/c^. These results 
give the best measurement of the y(4360) and y(4660). 

There is a faint evidence of y(4260) ^ tt+7i^^{2S), with i3(y(4260) 7r+7r-?/'(25)) • 
Fe+e- < 4.3 eV/c2 at the 90% confidence level, or S(y(4260) ^ TT+7r'ij{2S)) ■ Fg+e- = 
7.41^;^ eV/c^, depending on the interference between the y(4260) and y(4360). These 
numbers should be cornpared with the couplings of the y(4260) to n^Ti^J/tp in Ref. [3| and 
to K+K-J/ip in Ref. [3]. 
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